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2 Conservation Laws with Diffusion and Dispersion 

1. Introduction 

This paper is concerned with the convergence of smooth solutions u — w^' to 
the initial value problem for the nonlinear dispersive equation 

dtu + d^f{u)=ed^(3{d^u)-5dlu, (x,t)GlRx(0,oo), (1.1) 

with the initial condition 

u{x,0)=Uq (x), X G M, (1-2) 

where the parameters £ > and 6 > will converge to zero and Wq' is an ap- 
proximation of a given initial condition -uq • IR- ~^ IR- The flux / = f{u) and 
the (degenerate) viscosity (3 = /3(A) are given smooth functions satisfying certain 
assumptions to be listed shortly. When e = 0, the equation (1.1) reduces to the 
generalized Korteweg-de Vries (KdV) equation, the original one corresponding to 
the special flux f(u) = v? jl. On the other hand, when 5 = 0, we recover a non- 
linear degenerate parabolic equation. The flrst case was extensively studied from 
mathematical and numerical standpoints; see for instance [13, 14, 6] and the ref- 
erences therein (see also Section 2 for some background). For 5 = and under 
suitable assumptions, this equation was treated in [22] as a simplifled model of the 
pseudo viscosity approximation proposed by von Neumann and Richtmyer [24] and 
then studied in [26, 25] for numerical purposes. 

In this paper we prove, under an assumption on the relative size of the parameters 
e and 5, that the sequence {vF'^^ converges to an entropy solution u of the flrst 
order hyperbolic conservation law 

dtu + d^f{u)=^ (1.3) 

as e and 5^0. Throughtout the paper, we assume the following conditions on / 
and /9: 
(Ai) There are two constants Ci > and m > 1 such that 

|/'(w)| <Ci(l + |wr-^) forallwGlR. 

(A2) The function /3 = /9(A) is non-decreasing and satisfles /9(0) = 0. 

Observe that, from (A2), we can deduce /9(A)A > for all A G IR. The following 
assumption will also be in use: 
(Bi) There exist constants C2, C3, A^ > 0, and r > 1 such that 

C2 lAI^'^ < /3(A)A < C3 lAI^'^ for all |A| > A. 

The function /?, therefore, is at least quadratic and could vanish on a bounded 
interval. Two other assumptions on the function /3 will be required for certain 
results below: 
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(B2) There are constants C4 > and A^ > such that 

/3(A)A > C4 \Xf for all |A| > N. 

(B3) There are constants C5 > and r > 1 such that 

/9(A)A > C5 |A|3^ for all A e M. 

Observe that both (B2) and (B3) are refinements to the lower bound estimate 
in (Bi), and that (B2) is actually a consequence of (B3). For our main theorems, 
we shall assume either (Bi), or (Bi) and (B3). The condition (B2) will be relevant 
only in the course of the discussion of the proofs. One typical example of a function 
(3 satisfying (A2) and (Bi) (and also (B2), (B3)) is given by 

P{X) = \Xf--^X (r>l). (1.4) 

Recall that, for the Cauchy problem for (1.3), existence and uniqueness of an 
entropy solution were first proved by Kruzkov [16] (Section 2 below). The vanishing 
viscosity limit, which corresponds to 5 = and £ ^ 0, was studied by several 
authors in the special case of a linear (and therefore non-degenerate) diffusion 
term /3(A) = A; this activity started with Oleinik's [23] and Kruzkov's [16] works. 
For more general viscosity coefficients and under the assumptions (A2)-(Bi), the 
convergence was proved in [22]. 

The vanishing dispersion limit {e = and 5-^0) has also been widely investi- 
gated, see Lax-Levermore [20], [7, 30, 31], and the survey paper by Lax [19]. The 
main source of difficulty lies in the highly oscillating behavior of the solutions u^' , 
which do not converge in a strong topology. 

The complete equation (1.1) arises, at least for the linear diffusion /3(A) = A 
(the so called Korteweg-de Vries-Burgers equation), as a model of the nonlinear 
propagation of dispersive and dissipative waves in many different physical systems. 
For related studies, let us refer to [1, 2, 3, 10, 11]. Travelling waves were studied in 
[2] and [9] and (partial) analytical results can be found in [15] and [6]. For a study 
of numerical approximations with similar features as (1.1), see [4]. 

The vanishing dispersive and diffusion limit (both e and 6 tend to zero) was 
first studied by Schonbek [27]. She assumed /3(A) = A and used the compensated 
compactness method introduced by Tartar [29]. In particular, she proved that 
the sequence {u^'^} converges to a weak solution to (1.3) (but not necessarily the 
unique entropy one) under the assumption that either d = 0(£^) for f{u) = w^/2, 
or 5 = 0(£^) for arbitrary subquadratic fiux-functions / (i.e. take to = 2 in (Ai)). 
Our inequalities on e and 6 requires that viscosity dominates dispersion, which 
is expected since {u^'^} are known not to converge to a weak solution of (1.3) 
(a fortiori to the entropy one) when dispersion effects are dominant. 

Here we consider the general nonlinear equation (1.1) and we show that, under 

5 — rn, r- 

the assumption (Bi) and if 5 < C e^-"^ [m < 3), the solutions -u^' converge to 
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the entropy weak solution of (1.3) (Theorem 4.1). Under the assumptions (Bi) and 

5 — m 

(B3) and when 6 < Ce""^^-'"'-^ {m < 5 — 1/r, r > 1), the same result holds true 
(Theorem 4.2). Furthermore, our Theorem 4.3 below improves upon Schonbek's 
result in [27] in the case that /3(A) = A and / is subquadratic. 

The main tool necessary to deal with these singular limits is the concept of 
entropy measure- valued solution to the equation (1.3). Measure- valued solutions 
were introduced by DiPerna [7] (see Section 2) to represent the weak-* limits of 
the solutions of (1.1), According to DiPerna's theory, the strong convergence of an 
approximate sequence follows once one obtains 
(i) uniform bounds in L°°(0, T; L'^(IR)) for all T > with q > m (where m is given 

by (Ai); 
(ii) weak consistency of the sequence with the entropy inequalities and the strong 

consistency with the initial data (Cf. the conditions (2.5) and (2.6) of Section 

2). 

Such a framework can be used to establish as well the convergence of difference 
schemes as shown in [5]. Observe however that, because of the oscillatory behav- 
ior of w^' , no maximum principle is available for the equation (1.1) and the L^ 
(uniform) estimates are more natural. Therefore we actually use the L'^ version of 
DiPerna's result derived by Szepessy in [28]. 

One basic estimate that does not require any restriction on the parameters e 
and d is the so-called energy estimate based on the simplest quadratic function -u^ 
(Lemma 3.1). The derivation of the necessary L'^' estimates with larger exponents 
{q > 2) is based on the nonlinearity of the viscosity coefficient /3 and requires 
suitable restriction on s and 6. 

The problem studied in this paper is motivated by the numerical method pro- 
posed by Von Neumann and Richtmyer [24] (see also [26]) and based on the pseudo- 
viscosity approximation /3(A) = |A|A. An important advantage of the nonlinear 
(quadratic) diffusion term /3(A) is that it adds sharply localized viscosity near shocks 
and a small quantity (possibly zero if /3 vanishes for A < 0) elsewhere. The pseudo- 
viscosity idea was proposed by Von Neumann to stabilize an unstable scheme such 
as the Lax-Wendroff scheme [21, 18, 8]: a minimal amount of numerical viscosity 
is added in order to prevent both the formation of unphysical (entropy violating) 
shocks and the generation of highly oscillatory approximate solutions. 

2. Preliminaries 

This section contains short background material on Young measures, entropy 
measure valued (mv) solutions (Section 2.1), and dispersive equations (Section 2.2). 

2.1 Young Measures and Measure- Valued Solutions 

Following Schonbek [27], we describe a representation theorem for Young mea- 
sures associated with a sequence of uniformly bounded functions of L*^ where 
q G (1,00) is fixed in the whole of this subsection. The corresponding setting 
in L°° was first established by Tartar [29]. 
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Lemma 2.1. Let {uj} he a uniformly hounded sequence in L°°(]R_|_; L^(IR)). Then 
there exists a suhsequence {uj'} and a weak-k measurahle, mapping «/ : IR x IR+ -^ 
Prob (M) taking its values in the spaces of non-negative measures with total mass 
one (prohahility measures) such that, for all functions g G C(IR) satisfying 

g{u) = 0{\u\'') asl-ul^oo (2.1) 

for some r G [0, q), the following limit representation holds 



lim // g{uji{x,t))(p{x,t) dxdt — // / (/(A) dz/(a;^t)(A) (/)(x, t) dxdt 

j'^^ J JmxM+ JJjrxJR+Jm 

(2.2) 

for all (j) e C^ {M X 11+) . 

The measure-valued function i^(x,t) is a Young measure associated with the se- 
quence {uj}. The following result reveals the connection between the structure of 
1/ and the strong convergence. 

Lemma 2.2. Suppose that u is a Young measure associated with a sequence {%} 
that is uniformly hounded in L°°(1R+; L9(IR)). For u G L°°(IR+; L'^i'IR)), the 
following statements are equivalent: 

(i) limj^oo Uj = u in L°°(IR-|-; L\^^ (K-)) for some r G [1, q); 

(ii) v^x^t) = Su{x,t) for almost every (x,t) G IR x 1R+. 

In (ii) above, the notation du{x,t) is used for the Dirac mass defined by 



< Su(x,t),9{-) > (p{x,t)dxdt= // g{u{x,t))(p{x,t)dxdt 

MxTR+ J Jm.x]R+ 

for all g G C(IR) satisfying (2.1) and all (p e C^{Mx IR+). Following DiPerna 
[7] and [28], we now define the measure- valued solutions to the first order Cauchy 
problem 

dtu + dxf{u)=0, (2.3) 

u{x,0) =uo{x), xeM. (2.4) 

Definition 2.1. Assume that f satisfies the growth condition (2.1) and uq G 
L^(IR) n L^(IR). A Young measure v associated with a sequence {uj\, which is 
assumed to he uniformly hounded in L°°(IR-(-; L'^(IR)), is called an entropy mv so- 
lution to the prohlem (2.3)-(2.4) if 

at(z/(.), |A - k\) + a,(z/(.), sgn (A - k){f{\) - f{k))) < (2.5) 

in the distrihutional sense for all k eM, and, for all interval / C IR, 

lim—/ {u(^x,t)A^-uo{x)\)dxdt = 0. (2.6) 

T^0+ 1 Jq Jj 



6 Conservation Laws with Diffusion and Dispersion 

Remark. A function u E L°°(IR+; L^(IR) fl L'^(]R)) is an entropy weak solution 
to (2.3)-(2.4) in the sense of Kruzkov [16] if and only if the Dirac measure du{-) is 
an entropy mv solution. In the case p = +00, existence and uniqueness of such 
solutions were proved in [16]. The following results of entropy mv solutions were 
proved in [28]: Theorem 2.3 states that entropy mv solutions are actually Kruzkov's 
solutions. Theorem 2.4 states that the problem has a unique L^ solution. 

Theorem 2.3. Assume that f satisfies (2.1) and uq G L^(Sl) fl L'^(IR). Suppose 
that V is an entropy mv solution to (2.3)-(2.4). Then there exists a function w G 
L°°(IR+; Li(lR) n ^^(IR)) such that 

^{x,t) = Sw{x,t) (2.7) 

for almost every (x, t) G IR x IR_|_. 

Theorem 2.4. Assume that f satisfies (2.1) anduo G L-'^(IR) nL^(IR). Then there 
exists a unique entropy solution u G L°°(IR-|-; L^(M?) fl L^(IR)) to (2.3)-(2.4) which, 
moreover, satisfies 

I|w(-,^)I|l'-(]R) < ||wo||l'-(]r) (2.8) 

for almost every t G 1R+ and for all r E [1, q] . Moreover the measure-valued mapping 
^{x,t) = ^u{x,t) is the unique entropy mv solution of the same problem. 

Combining Theorems 2.3 and 2.4 and Lemma 2.2 we obtain our main convergence 
tool. 

Corollary 2.5. Assume that f satisfies (2.1) and uq G L-'^(IR) fl L'^(IR). Let {uj} 
he a sequence of functions that are uniformly hounded in L°°(1R_|_; L'3'(1R)) for q>l, 
and let f he a Young measure associated with this sequence. If f is an entropy mv 
solution to (2.3)-(2.4), then 

lim Uj - w in L°°(1R+; L[„^(IR)) (2.9) 

for all r G [1,?), where u G L°°(1R+; L'^(IR)) is the unique entropy solution to 
(2.3)-(2.4). 

2.2 Nonlinear Dispersive Equations. 

Consider the fully nonlinear, KdV-type equation in one space variable 

dtu + F{u, d^u, dlu, dlu) = (2.10) 

for (x,t) G IR X (0, T) together with the Cauchy data 

u{x,0)=uo{x) for a; G IR. (2-11) 

In recent years this class of problems has been extensively studied; see for in- 
stance [13, 14, 15] and the references cited therein. In particular, under suitable 
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assumptions, these equations enjoy a gain of regularity of the solutions with respect 
to their initial data (Cf. [6, 12, 17] and the references above). We observe however 
that a complete theory of global existence remains to be developped. Presenting a 
complete review is far beyond the aim of this presentation. Here we need only recall 
a result -that is sufficient to deal with (1.1)- of local existence and uniqueness ot 
the smooth solutions to the equation (2.10). More details and proofs can be found 
in the paper by Craig-Kappeler-Strauss [6] (see also [15]). 

Using the notation U = {uq, ui,U2, ws), the assumptions on the smooth function 
F = F{U) are as follows: 
(Hi) There exists d > such that 



du^F{U) > 5 > for all t/ G M 
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H2) du^F{U) < for all U e IR^. 

The equation (1.1) does satisfy these hypotheses provided d > and (A2) holds. 

Theorem 2.6 (Uniqueness). Let T > be fixed and assume F satisfies (Hi)-(H2). 
For any uq G if^(IR) there is at most one solution u G L°°((0, T); i7'^(IR)) to the 
problem (2.10)-(2.11). 

Theorem 2.7 (Existence). Assume F satisfies (Hi)-(H2). Let N be an integer 
> 7 and let Cq > be a given constant. There exists a time T > 0, depending only 
on Co, such that for all uq G H^(]R) with ||'Uo||h'' < Cq, there exists at least one 
solution u G L°°((0,T); if^(lR)) to the problem (2.10)-(2.11). 

Observe that the setting above based on a Sobolev norm of relatively high order 
(ilf^(IR)) is the optimal result provided by the current techniques of analysis of 
dispersive equations. The global existence of smooth solutions to (1.1)-(1.2) appears 
to be an open problem. In the following we tacitly restrict attention to a time 
T* G (0, cxd] chosen such that the problem (1.1)-(1.2) is well-posed in the strip 
Mx (0,T,). 

3. A Priori Estimates 

In this section we consider a sequence {u^'^} of smooth solutions of (1.1)-(1.2) 
that vanish at infinity. We assume also that the initial data {wq' } are smooth 
functions with compact support, and are uniformly bounded in L^(IR) fl L'^(IR) for 
a suitable q > 1. In what follows, whenever it does not lead to confusion, we omit 
the indices e and S. Similarly all constants are denoted by C, Ci, ■ ■ ■ . We begin 
with the natural energy estimate based the quadratic function u'^/2. (Arbitrary 
convex functions could not be easily used here because of the dispersive term.) 

Lemma 3.1. For any T > 0, it holds that 

u^{x,T) dx + 2e / / /3{ux{x,t))uj;{x,t) dxdt = / UQ{x)dx. (3.1) 

IR Jo Jm Jm. 
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Proof. We multiply (1.1) by u and integrate in space. Integrating by parts we 
obtain 



A 1 

dt 2 



/ u^ d. + / „/'(„)«. d. = -e / ^Mu^ i.-5[ ««... d.. (3.2) 
Jtr Jtr Jtr Jtr 



The second terms on both sides of (3.2) vanish identically, since we can write these 
terms in a conservative form: 

uf'{u)un, = {G{u))^ 

with G' = uf'{u) and 

The estimate (3.1) follows therefore from (3.2). D 

At this stage, in view of the assumption (Ai) and Lemma 3.1, we have a uniform 
control of w in L°°(IR+; L2(IR)) and e I3{u^)u^ in Li((0,T) x IR) for every T > 0. 
Our next aim is to derive a (non-uniform) estimate of u in L°° norm, which will be 
uniform in e but not in 5. We first provide a second energy- type estimate. 

Lemma 3.2. Let F he defined by F'{u) = f{u). For every T > 0, we have 



- / ul{x,T)dx- / F{u{x,T))dx + e6 / / ul^l3'{u:,)dxdt 
^ Jm Jtr Jo Jtr 

= o / ^o,x(^)dic- / F{uo{x))dx + e / /'(w)/3(w^)w^ dxdt. 

^ Jm ' Jm Jo Jtr 

Proof. Multiplying (1.1) by f{u) + Suxx we obtain the equality 

0= / {F{u)t + H{u)x - ef{u)(3{ux)x + Sf{u)uxxx} dx 
Jm. 



(3.3) 



IR 

2 o'/ 



+ / {SUxxUt + SUxxf{u)x- Su^^l3'{Ux) +SUxxxUxx}dx 

Jm. 
where we have set H' = ff. After integration by parts in space, we have 

/F{u) dx — -u^ dx + / ef'{u)/3{ux)ux dx — e6 / u^^/3'{ux)dx = 0, 
2 Jm Jtr 



d 
dt 



which gives the desired conclusion. D 

Using the bound for v5 Ux in L°°(IR+; L^(IR)) that follows from Lemma 3.2, we 
are now able to estimate u in the L°° norm. 
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Lemma 3.3. If m < 5 in the assumption (Ai), there exists a constant C > such 
that 

\u{x,t)\<C6~^^ (3.4) 

for all {x,t) eMx 1R+. 

Proof. In view of (A2) and for all w G M, we have 

\F{u)\ <C{l + \ur+^). 

The main idea is to use (3.3) to control du'^ in terms of F{u), the latter being 
estimated from the above growth condition. We deduce from (3.3) that 

- / u1{x,T)dx + e6 / u'?^^ /3' (ux) dx dt 

2 J]R Jo J-R 

= o[ «o,x(a^) dx+ f {F{u{x, T)) - F{uQix))) dx + e f f f'{u)(3Mux dx dt 
^ Jtr ' Jm Jo J 

<C + C||^-,T)||--i^J ||^.,T)||i.(n,)+e^^ j||^/3(tx.)z.. dxdt 

where we used (A2) as well. (We do not explicitly write the terms involving the 
initial data as it is assumed to be a smooth function.) Therefore in view of (3.1), 
we arrive at 



2 



f ulix,T)dx + e5 f f ul,(3'Mdxdt<C +C\\u{-,T)\\'J^-I^y 

JjR Jo Jm 

Using the Cauchy-Schwartz inequality, we obtain 



'u(x,t)p<2/ \u{y,t)ux{y,t)\dy 



2 

V6' 



< ^^ll«(-,^)l|L2(]R)V^||Wx(-,t)||L2(IR) 



_C „ /. . ^_, \l/2 

Hence, for alH > 0, 



<^ll«0||L2(]R)(l+||tx(-,t)||^,(i))" 



\\<MUm ^ f (1 + ll^(-'^)lir^?]R)) • (3-5) 

Now, since to < 5, the growth of the left hand side of (3.5) exceeds the growth 
of the right hand side. So we can check that (3.5) implies 

lk(-,^)l|L-(]R) < max ( 1, ( -7- 
< C 5 5-™ . 
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Namely, setting y = ||w(-, t)||Loo, we have for y > 

y' < J {I + y"^-') . (3.6) 

If y < 1 we have the conclusion. Otherwise, suppose we had y > (^) ^~"" , then 
we would deduce that 

y'>^-^y'--'>j{i + y'-'), 

which would contradict the inequality (3.6). D 

In view of the proof of Lemma 3.3 we also state the following result: 
Lemma 3.4. For any T > we have 

-T 



-I ul{x,T)dx + e / ul^/3'{u^)dx dt<CS~^. (3.7) 

2 Jm. Jo Jm 



Finally we derive uniform bounds in L°°((0, T); L'^(IR)) with q < 5. The follow- 
ing result provides a uniform estimate in the Lebesgue space L^(IR) (so improving 
upon the L^ bound in Lemma 3.1) by taking advantage of the nonlinearity property 
of (3 (for large values of A = Ux) as stated in (B2). 

Proposition 3.5. Assume that the assumption (B2) holds and m < 3 in the as- 
sumption (Ai). There exists a constant C > 0, which depends only on the initial 
data, such that, for all d and e small enough and for every T > 0, we have 

sup lk(-,t)||ia(n,) < C (T + e-'6l^) . (3.8) 



Proof. Set r]{u) = \u\^. Multiply (1.1) by r]'{u) and integrate on IR x (0, T). It 
follows easily that 

ri{u{x,T)) dx + e j j ri"{u)/3{uj;)uj; dxdt 

(3.9) 
ri{uQ{x))dx / / rj'" {u)u\ dx dt. 



IR ^0 ^IR 



IR 2 Jo JIR 



On the other hand, according to (B2), there exists a constant C such that for all 
A GIR 

|Ap<C(l + /3(A)A). 
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Using the latter in the energy estimate (3.1), we are able to control the second term 
in the right hand side of (3.9): 



V"'{u)ul dxdt 



<Ci / \u\'^\u^fdxdt 

Jo Jm. 

<C2 / / \u\^ {l + /3{u^)u^) dx dt 
Jo Jm 

C 
<C3T sup ||w(t)||i,2(n^) + -||'u|||^(ij^x(o^^)) 



te(o,T) 



e 



<CT +C il + e-^d^ 



where the latter inequality follows from the L°" bound in Lemma 3.3. Returning 
to (3.9), we obtain (3.8). D 

Under the stronger assumption (B3), we have a sharper estimate. 

Proposition 3.6. Assume that (B3) holds for a given r > 1 and that m < q = 
5 — 1/r in (Ai). There exists a constant C > 0, which depends only on the initial 
data, such that, for all d and e small enough and every T > 0, 



sup \H;t)\\U^)<C (T'-'/^ + s-'/^d-i^). 
te{o,T) '' \ J 



(3.10) 



Proof. As in the proof of Proposition 3.5, we consider the formula (3.9) but now 
with r\{u) = \u\^ with g = 5 — 1/r. The assumption (B3) yields, using (3.1) and 
(3.4), 



ri"'{u)u^ dxdt 



^IR 



<C I I \u\'^-^\u^\^ dx dt 
'0 Jtr 

<c \ I [ i^r'^^-^Mxdt 

'0 Jm 



1/r' 



1/r 



l3{ux)Ux dxdt 



JWi 



<Ce / |hllLoo(iR)7^ ^ sup ||w(-,t)||/2(iR) 



te(o,T) 



< CT^/'^'e-^/'^^" 



r(5 — m) 



where 1/r + 1/r' = 1 and since r'{q — 3) = 2 + r' /r. Now (3.10) follows from the 
formula (3.9). D 

Observe that the estimates (3.8) and (3.10) hold only thanks to the nonlinear 
form of the viscosity term. Note that when 5 = 0, these estimates reduce to bound 
that are uniform in e, which is expected for conservation laws with viscosity but no 
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dispersion. On the other hand, most of our estimates blow-up when taking e = 
and provide no control of L*^ norms. Observe however that the L^ bound in Lemma 
3.1 is uniform for 6 arbitrary and £ ^ 0. Loosely speaking, this is also the case 
of the estimate in Proposition 3.6 if r could be chosen to be r = oo, in which case 
(3.10) becomes a uniform L^ estimate. Such a choice of r is not allowed however, 
cf. (Bs). 

For the sake of completeness, we finally state an analogous estimate for linear 
diffusions which was proved by Schonbek [27]. 

Proposition 3.7. Let /3(A) = A and take m = 2 in (Ai). For any T > 0, there 
exists a constant Ct > 0, which depends only on the initial data, such that for 

sup ||'u(-,t)||L4(iR) < Ct- (3.11) 

t6(0,T) 



4. Convergence Results 

In Section 3, we have established several uniform bounds for the sequence {u'^'^} 
of solutions to the Cauchy problem (1.1)-(1.2) under certain assumptions on the 
functions /, /?, and the parameters e and 5. Assume again that the initial data 
Uq are smooth with compact support and that there exists a limiting function 
uq G L^(IR) n L'^(IR) and a suitable q> 1 (specified below) such that, if 5 = 0(e), 

\iuvu/ = uq inLi(IR)nL5(lR). (4.1) 

Returning to the proofs of Section 3, it is not hard to see that the following con- 
ditions on the initial data are sufficient for the estimates therein to hold uniformly 
with respect to a class of initial data: 

II Ei^ll I II s,5\\ rl/2 II e,5 \\ ^ r^ 

IFo IU2(]R) + \\Uq ||l9(]R)0 ' ||Wo,xlU2(]R) < C 

In this section we prove the strong convergence of the sequence w^'^. 

Theorem 4.1. Assume that (Bi) holds and m < 3 in (Ai). Letu^'^ be a sequence 
of smooth solutions to (1.1)-(1.2) on IRx (0, T) (for a given T > 0), which vanish at 
infinity and are associated with initial data satisfying (4.1) with q = 5. If there is a 

5 — m r- 

constant C > such that S < C e^-"^ , then the (whole) sequence u^' converges to a 
function u G L°°((0, T); L^(IR)), which is the unique entropy solution to (2.3)-(2.4). 

Theorem 4.2. Assume that (Bi) and (B3) hold and to < 5 — 1/r in (Ai). Let 
u^'^ be a sequence of smooth solutions to (1.1)-(1.2) on IR x (0, T) (for a given 
T > 0), which vanish at infinity and are associated with initial data satisfying (4.1) 

5 — m 

with q = b — 1/r. If there is a constant C > such that 6 < C£^(5-m)-i ^ then 
the (whole) sequence converges to a function u G L°°((0, T); L^(IR)), g = 5 — 1/r, 
which is the unique entropy solution to (2.3)-(2.4). 

Let us give an analogous statement in the case /9(A) = A, which improve upon 
[27] (Cf. Theorem 5.1 therein). 
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Theorem 4.3. Assume m = 2 in (Ai) and let (3{X) = X. If d < C e^ , the whole se- 
quenceu^'^ of solutions o/(l.l)-(1.2) converges to afunctionu G L°°((0, T); L'^(]R)), 
which is the unique entropy solution to (2.3)-(2.4). 

Similar convergence results can be proven for the case to = 3 or for some special 
flux- functions along the lines of what was done in [27] (Cf. Sections 4 and 5 therein). 
Theorems 4.1-4.2 follow easily by simply using the following general result and the 
L'^ bounds derived in Section 3, Proposition 3.5 and Proposition 3.6 respectively. 

Theorem 4.4. Assume that (Bi) holds. Letu^'^ be a sequence of smooth solutions 
to (1.1)-(1.2) onJRx (0, T) (for a given T > 0) associated with initial data satisfying 
(4.1). If the sequence is uniformly bounded in L°°((0, T); L'3'(IR)) forq > m and d = 
o{e^/'^), then the (whole) sequence converges to a function u G L°°((0, T); L'^(IR)), 
which is the unique entropy solution to (2.3)-(2.4). 

Proof of Theorem 4 ■4- First of all let us establish that, for any convex function 
rj = rj{u) such that rj' , r]" , rj'" are uniformly bounded on IR, we have 

A^'^ ^dt'n{u^'^)+d^Q{u^'^)^0 mV'{MxM+), (4.2) 

where Q' = f'r( . To begin with, observe that 

A'^'^ = e{rl{u)^{u^))^ - er]"{u)(3{u^)u^ - 5{r)'{u)u^^)^ + 5r}"{u)u^u^^ 

= Ti+T2 + T3+T4. 

To estimate Ti, we use the assumption (Bi), which implies |/3(A)| < C (l -|- 
|^|3r-i^ for aU A. For any given ^ G C^(IR x (0, T)), ^ > 0, and for p = ^ > 1 
{p' being the conjugate exponent of p) , we get 



i/p 



(Ti,^)= II ee,T]'{u)(3{u,)dxdt 

< C'^||^x||l1(]Rx(0,T)) +C'£||^x||lp'(]Rx(0,T)) ( 

SO using the derivative estimate in (3.1) and (Bi) again: 

(Ti,6') <C(£ + £i/3'^) < Ce^/3^ (4.3) 

The second term T2 is nonpositive, namely 

{T2,e) = - ffe7]"iu)/3{u^)uj dxdt < 0. (4.4) 
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To estimate T3, we use the energy estimate (3.1) and the fact that (3 is at least 
quadratic. We write 



(T3, 6') = (5 / / Oj:!]' {u)uxj: dx dt 

'lRx(0,T) 

^ 11 Ox{{'n'{u)u^)^-r]"{u)ul) dxdt 

'lRx(0,T) 



< -S Oj:^r]'{u)uj: + S \9^\uldxdt 

J JtRx(0,T) J Jmx{0,T) 

<CS+ C6 ( ff \u:,f'^dxdt] + Ce-2/3'^, 

\J J supp 9 J 

where supp6' denotes the support of the function in IR x (0, T), thus 
(T3, 6*) < C (5 (1 + C e-'^l^'' + £-^/3") < C (1 + C e-'^l^'' 
Finally we deal with T4 as follows: 



(4.5) 



1(^4,^)1 



S 1 1 rj '{u)ux UxxO dx dt 
6 



1 



9W'{n)ui + d II 9v"'{u)-ut 



<Cd 



\ux\ + \uj:\ dxdt 



supp 



SO 



\{T^,e)\<C6 ('e-2/3'^ + £-1/'^) < C e-^l' 



(4.6) 



Therefore, if 5 = o{e^''^\ (4.2) follows immediately from the estimate (4.3)-(4.6). 
To apply Corollary 2.5 we have to show that (2.5) and (2.6) are satisfied for a Young 
measure v associated with the sequence u'^'^ . It is a standard matter to deduce, for 
all convex entropy pairs. 



at(z/(.), 77(A)) + a,(z/(.),Q(A))<o 

from the convergence property (4.2). The inequality (2.5) for all /c G IR then follows 
by using a standard regularization of the function \u — k\. Concerning the initial 
data and in order 

to establish (2.6), we now combine the entropy inequalities and the weak consis- 
tency property as was suggested by DiPerna [7] . We follow the detailled arguments 
given in [28]. 
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Consider the function g{X) = |A|^ for 1 < r < min (2,q), and set 

G(A, Ao) = g{X) - g{Xo) - <7'(Ao)(A - Aq) 

r(r-l) (A-Ao)2 (4.7) 



> 



2 (i + |A| + |Ao|)2-- 



Let / C IR be be a closed and bounded interval. Using the Jensen inequality, the 
Cauchy-Schwartz inequality, and the above convexity inequality (4.7), it is easily 
checked that 



1 "^ 



T.,0 J I 



{^{x,t), |A-'Uo(a:^)|)dxdt 

T . nV2 (4.8) 

T 



<Ci \— I I (z/(^,t), (^(A, uo{x))) dx dt 



Let {i^n} G Cq°(SI) be a sequence of test-functions such that 

lim i/jn = g'{uo) = in L'' (R), 

where 1 = 1/r + 1/r' . Using the uniform bound in L*^ available for the sequence 
{uq^}: we get 

/ {iy{x,t),G{X,uo{x)))dxdt 
Ji 

T 



(4.9) 

. I |Uo| ux' 

+ '2T \\uo\\L-{]R)\\g' (Uo) -i^n\\L-'{]R)- 



- / {'^{x,t),uo- X)tlJn dxdt + T \uo\'^dx 

Jo Jm Jm-^i 



Taking an increasing sequence of compact sets Ki covering IR, i.e. such that / C 
Ki C K2 C . . . and Ui^i ^i = ^^ we have 

{^{x,t),G{\uo{x)))dxdt < / (z/(^^t),G(A,'Uo(x)))dxdt, 

Jl Jo JKi 

which, together with (4.9) where / is replaced by Ki, yields 



1 '■^ 



1 



Ji 

T 



{^{x,t),Q{>^,uo{x)))dxdt 



<— / / {iy{x,t),Uo{x)- X)'l{jn dxdt + 2\\uo\\L-{M)\\9'{uo)-i^n\\L-'(M)^ 

-^ '0 Jm 



(4.10) 
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Since 

|2 



lim / \uo\ dx — 0. 

Therefore, in view of (4.8) and (4.10), the strong consistency property (2.6) will 
be established if we show that 



lim — / / (z/(^^t),wo(a^) - A)V'n dx dt < (4.11) 



tZo+ T 
for all n G N. By definition of the Young measure (Cf. (2.2)), we have 

7f / Wix,t),uo{x)-X)^n da;dt= lim — / / (uo{x) -u^'^{x,t)) i/jn dxdt. 

J- Jo Jm ^'^^^ J- Jo J 

On the other hand, we can write 

= / [uo{x) -UQ^{x,t))iJri{x)dx - — / / ( / dsu'''^{x,s)ds\^n{x)dxdt 

Jm. ^ Jo Jm \Jo J 

= A + B. 

The term A tends to zero as £ -^ in view of the weak consistency property (4.1). 
Furthermore, by arguing as in the derivation of (4.3), we have 

^--fl J(l {-d.f{u''')+ed,{P{ul''))-6dlu'Ai;^{x)dxdt 

= -^[ [ [ {f{u'^')d,ijr.-e(3{uf)d,ijr. + 5u'^'dlijr.)dsdxdt 

^ Jo JtrJo 

<CnT. 

This leads to the inequality (4.11). The proof of Theorem 4.4 is completed. D 

The proof of Theorem 4.3 is based on slightly modified estimates in the inequal- 
ities (4.5) and (4.6), which can be derived by arguing as in [27]. The details of the 
proof are omitted. 
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